Genomic selection for morphological and yield-related traits using genome-wide SNPs in oil palm.

Oil palm is the most important oil crop worldwide. Colombia is the fourth largest producer, primarily relying on production from interspecific hybrids, derived from crosses between Elaeis oleifera and Elaeis guineensis (OxG). However, conventional breeding can take up to 20 years to generate a new variety. Therefore, reducing the breeding cycle while improving the genetic gain for complex traits is desirable. Genomic selection (GS) is an approach with the potential to achieve this goal. In this study, we evaluated 431 F1 interspecific hybrids (OxG) and 444 backcrosses (BC1) for morphological and yield-related traits. Genomic predictions were performed with the G-BLUP model using three different population datasets for training the model: the same population (TRN1), the other population (TRN2), and both populations (TRN1+2). Higher multi-family prediction accuracies were obtained for foliar area (0.3 in OxG) and trunk height (0.47 in BC1) when the model was trained with TRN1. Single-family prediction accuracies were lower in the OxG compared to BC1 families for traits such as trunk diameter, trunk height, bunch number, and yield using TRN1. Conversely, lower prediction accuracies were obtained for most traits when the model was trained using TRN2 (< 0.1). Multi-trait models showed a substantial increase of the predictions for traits such as yield (0.22 for OxG and 0.44 for BC1), because of the genetic correlations between traits. The results herein highlighted the potential of GS for parental selection in OxG and BC1 populations, but further studies are required to improve the models to select individuals by their genetic value. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-022-01341-5.

Comparative transcriptomic analysis reveals key components controlling spathe color in Anthurium andraeanum (Hort.).

Anthurium andraeanum (Hort.) is an important ornamental in the tropical cut-flower industry. However, there is currently insufficient information to establish a clear connection between the genetic model(s) proposed and the putative genes involved in the differentiation between colors. In this study, 18 cDNA libraries related to the spathe color and developmental stages of A. andraeanum were characterized by transcriptome sequencing (RNA-seq). For the de novo transcriptome, a total of 114,334,082 primary sequence reads were obtained from the Illumina sequencer and were assembled into 151,652 unigenes. Approximately 58,476 transcripts were generated and used for comparative transcriptome analysis between three cultivars that differ in spathe color ('Sasha' (white), 'Honduras' (red), and 'Rapido' (purple)). A large number of differentially expressed genes (8,324), potentially involved in multiple biological and metabolic pathways, were identified, including genes in the flavonoid and anthocyanin biosynthetic pathways. Our results showed that the chalcone isomerase (CHI) gene presented the strongest evidence for an association with differences in color and the highest correlation with other key genes (flavanone 3-hydroxylase (F3H), flavonoid 3'5' hydroxylase (F3'5'H)/ flavonoid 3'-hydroxylase (F3'H), and leucoanthocyanidin dioxygenase (LDOX)) in the anthocyanin pathway. We also identified a differentially expressed cytochrome P450 gene in the late developmental stage of the purple spathe that appeared to determine the difference between the red- and purple-colored spathes. Furthermore, transcription factors related to putative MYB-domain protein that may control anthocyanin pathway were identified through a weighted gene co-expression network analysis (WGCNA). The results provided basic sequence information for future research on spathe color, which have important implications for this ornamental breeding strategies.

Combining transcriptome analysis and GWAS for identification and validation of marker genes in the Physalis peruviana-Fusarium oxysporum pathosystem.

Vascular wilt, caused by the pathogen Fusarium oxysporum f. sp. physali (Foph), is a major disease of cape gooseberry (Physalis peruviana L.) in Andean countries. Despite the economic losses caused by this disease, there are few studies related to molecular mechanisms in the P. peruviana-Foph pathosystem as a useful tool for crop improvement. This study evaluates eight candidate genes associated with this pathosystem, using real-time quantitative PCR (RT-qPCR). The genes were identified and selected from 1,653 differentially expressed genes (DEGs) derived from RNA-Seq analysis and from a previous genome-wide association study (GWAS) of this plant-pathogen interaction. Based on the RT-qPCR analysis, the tubuline (TUB) reference gene was selected for its highly stable expression in cape gooseberry. The RT-qPCR validation of the candidate genes revealed the biological variation in their expression according to their known biological function. Three genes related to the first line of resistance/defense responses were highly expressed earlier during infection in a susceptible genotype, while three others were overexpressed later, mostly in the tolerant genotype. These genes are mainly involved in signaling pathways after pathogen recognition, mediated by hormones such as ethylene and salicylic acid. This study provided the first insight to uncover the molecular mechanism from the P. peruviana-Foph pathosystem. The genes validated here have important implications in the disease progress and allow a better understanding of the defense response in cape gooseberry at the molecular level. Derived molecular markers from these genes could facilitate the identification of tolerant/susceptible genotypes for use in breeding schemes.

Optimization of the genotyping-by-sequencing SNP calling for diversity analysis in cape gooseberry (Physalis peruviana L.) and related taxa.

A robust Genotyping-By-Sequencing (GBS) pipeline platform was examined to provide accurate discovery of Single Nucleotide Polymorphisms (SNPs) in a cape gooseberry (Physalis peruviana L.) and related taxa germplasm collection. A total of 176 accessions representing, wild, weedy, and commercial cultivars as well as related taxa from the Colombian germplasm bank and other world repositories were screened using GBS. The pipeline parameters mnLCov of 0.5 and a mnScov of 0.7, tomato and potato genomes, and cape gooseberry transcriptome for read alignments, were selected to better assess diversity and population structure in cape gooseberry and related taxa. A total of 7,425 SNPs, derived from P. peruviana common tags (unique 64 bp sequences shared between selected species), were used. Within P. peruviana, five subpopulations with a high genetic diversity and allele fixation (HE: 0.35 to 0.36 and FIS: -0.11 to -0.01, respectively) were detected. Conversely, low genetic differentiation (FST: 0.01 to 0.05) was also observed, indicating a high gene flow among subpopulations. These results contribute to the establishment of adequate conservation and breeding strategies for Cape gooseberry and closely related Physalis species.

Genome-Wide Association Study Reveals Novel Candidate Genes Associated with Productivity and Disease Resistance to Moniliophthora spp. in Cacao (Theobroma cacao L.).

Cacao (Theobroma cacao L.), the source of chocolate, is one of the most important commodity products worldwide that helps improve the economic livelihood of farmers. Diseases like frosty pod rot caused by Moniliophthora roreri and witches' broom caused by Moniliophthora perniciosa limit the cacao productivity, this can be solved by using resistant varieties. In the current study, we sequenced 229 cacao accessions using genotyping-by-sequencing to examine the genetic diversity and population structure employing 9,003 and 8,131 single nucleotide polymorphisms recovered by mapping against two cacao genomes (Criollo B97-61/B2 v2 and Matina 1-6 v1.1). In the phenotypic evaluation, three promising accessions for productivity and 10 with good tolerance to the frosty pod rot and witches' broom diseases were found. A genome-wide association study was performed on 102 accessions, discovering two genes associated with productivity and seven to disease resistance. The results enriched the knowledge of the genetic regions associated with important cacao traits that can have significant implications for conservation and breeding strategies like marker-assisted selection.

Genome-wide association study (GWAS) for morphological and yield-related traits in an oil palm hybrid (Elaeis oleifera x Elaeis guineensis) population.

BACKGROUND: The genus Elaeis has two species of economic importance for the oil palm agroindustry: Elaeis oleifera (O), native to the Americas, and Elaeis guineensis (G), native to Africa. This work provides to our knowledge, the first association mapping study in an interspecific OxG oil palm population, which shows tolerance to pests and diseases, high oil quality, and acceptable fruit bunch production. RESULTS: Using genotyping-by-sequencing (GBS), we identified a total of 3776 single nucleotide polymorphisms (SNPs) that were used to perform a genome-wide association analysis (GWAS) in 378 OxG hybrid population for 10 agronomic traits. Twelve genomic regions (SNPs) were located near candidate genes implicated in multiple functional categories, such as tissue growth, cellular trafficking, and physiological processes. CONCLUSIONS: We provide new insights on genomic regions that mapped on candidate genes involved in plant architecture and yield. These potential candidate genes need to be confirmed for future targeted functional analyses. Associated markers to the traits of interest may be valuable resources for the development of marker-assisted selection in oil palm breeding.

Virus-induced gene silencing (VIGS) in cape gooseberry (Physalis peruviana L., Solanaceae) / Silenciamento de genes induzidos por vírus (VIGS) em physalis (Physalis peruviana L., Solanaceae) / Silenciamiento génico inducido por virus (VIGS) en uchuva (Physalis peruviana L., Solanaceae)

Abstract Cape gooseberry (Physalis peruviana, L.) is a herbaceous plant belonging to the Solanaceae family that produces an edible berry appreciated for its nutraceutical and pharmaceutical properties. Its production is often limited by diseases and reproducible fruit quality. Recent studies have reported genes associated with fruit quality and resistance response to the root-infecting fungus Fusarium oxysporum f. sp. physali (Foph,) which causes vascular wilt. In order to standardize a method to validate the biological function of candidate genes in the non-model species P. peruviana, we tested the robust approach in reverse genetics, virus-induced gene silencing (VIGS). In this study, we validated and optimized VIGS using an insert of the phytoene desaturase (PDS) gene in a silencing viral vector generated from tobacco rattle virus (TRV). Leaves infiltrated with Agrobacterium (GV3101 strain) showed photo-bleached segments, which were distinctive for PDS suppression at 7 days post-infection (dpi). More than half of the treated plants showed photo-bleaching, indicating an efficiency rate of 50 % of the VIGS protocol. The results of this study showed that VIGS can be used for future functional gene characterization implicated in the immune response, disease resistance and fruit quality in capegooseberry.

Resumo A physalis (Physalis peruviana, L.) é uma planta herbácea pertencente à família Solanaceae, que produz uma baga comestível apreciada por suas propriedades nutracêuticas e farmacêuticas. Sua produção com frequência se vê limitada devido a enfermidades e baixa reprodutibilidade na qualidade do fruto. Estudos recentes reportaram genes associados com a qualidade do fruto e com a resposta de resistência ao fungo radicular Fusarium oxysporum f. sp. physali (Foph.), que causa esmorecimento vascular. Com a finalidade de padronizar um método para validar a função biológica de genes candidatos na espécie não-modelo P. p ruviana, avaliamos uma aproximação robusta em genética invertida, o sil nciamento de genes induzidos por vírus (VIGS). Neste estudo, validamos e otimizamos o VIGS usando um inserto da fitoeno desaturase (PDS) em um vetor viral de silenciamento produzido a partir do vírus do chocalho do tabaco (TRV). As folhas infiltradas com Agrobacterium (cepa GV3101) mostraram segmentos fotobranqueados, que foram distintivos para a supressão de PDS a 7 dias pós-infecção (dpi). Mais da metade das plantas tratadas mostraram fotobranqueamento, o que indica uma taxa de eficiência de 50 % do procotolo VIGS. Os resultados de este estudo mostraram que o VIGS pode ser usado em caracterizações futuras de genes funcionais implicados na resposta imune, na resistência a enfermidades e na qualidade do fruto de physalis.

Resumen La uchuva (Physalis peruviana, L.) es una planta herbácea perteneciente a la familia de las solanáceas, que produce una baya comestible apreciada por sus propiedades nutracéuticas y farmacéuticas. Su producción con frecuencia se ve limitada debido a enfermedades y a falta de reproducibilidad en la calidad del fruto. Estudios recientes han reportado genes asociados con la calidad del fruto y con la respuesta de resistencia al hongo radicular Fusarium oxysporum f. sp. physali (Foph,), que causa marchitamiento vascular. Con el fin de estandarizar un método para validar la función biológica de genes candidatos en la especie no-modelo P. peruviana, evaluamos la aproximación robusta en genética inversa, el silenciamiento génico inducido por virus (VIGS). En este estudio, validamos y optimizamos el VIGS usando un inserto de la fitoeno desaturasa (PDS) en un vector viral de silenciamiento producido a partir del virus del cascabeleo del tabaco (TRV). Las hojas infiltradas con Agrobacterium (cepa GV3101) mostraron segmentos fotoblanqueados, que fueron distintivos para la supresión de PDS a 7 días pos-infeccion (dpi). Más de la mitad de las plantas tratadas mostraron fotoblanqueo, lo cual indica una tasa de eficiencia del 50 % del protocolo VIGS. Los resultados de este estudio mostraron que el VIGS se puede usar en futuras caracterizaciones de genes funcionales implicados en la respuesta inmune, la resistencia a enfermedad y la calidad del fruto en la uchuva.

Ancestry informative alleles captured with reduced representation library sequencing in Theobroma cacao.

As the source of chocolate, cacao has become one of the most important crops in the world. The identification of molecular markers to understand the demographic history, genetic diversity and population structure plays a pivotal role in cacao breeding programs. Here, we report the use of a modified genotyping-by-sequencing (GBS) approach for large-scale single nucleotide polymorphism (SNP) discovery and allele ancestry mapping. We identified 12,357 bi-allelic SNPs after filtering, of which, 7,009 variants were ancestry informative. The GBS approach proved to be rapid, cost-effective, and highly informative for ancestry assignment in this species.

Association Study Reveals Novel Genes Related to Yield and Quality of Fruit in Cape Gooseberry (Physalis peruviana L.).

Association mapping has been proposed as an efficient approach to assist plant breeding programs to investigate the genetic basis of agronomic traits. In this study, we evaluated 18 traits related to yield, (FWP, NF, FWI, and FWII), fruit size-shape (FP, FA, MW, WMH, MH, HMW, DI, FSI, FSII, OVO, OBO), and fruit quality (FIR, CF, and SST), in a diverse collection of 100 accessions of Physalis peruviana including wild, landraces, and anther culture derived lines. We identified seven accessions with suitable traits: fruit weight per plant (FWP) > 7,000 g/plant and cracked fruits (CF) < 4%, to be used as parents in cape gooseberry breeding program. In addition, the accessions were also characterized using Genotyping By Sequencing (GBS). We discovered 27,982 and 36,142 informative SNP markers based on the alignment against the two cape gooseberry references transcriptomes. Besides, 30,344 SNPs were identified based on alignment to the tomato reference genome. Genetic structure analysis showed that the population could be divided into two or three sub-groups, corresponding to landraces-anther culture and wild accessions for K = 2 and wild, landraces, and anther culture plants for K = 3. Association analysis was carried out using a Mixed Linear Model (MLM) and 34 SNP markers were significantly associated. These results reveal the basis of the genetic control of important agronomic traits and may facilitate marker-based breeding in P. peruviana.

Colombia a Source of Cacao Genetic Diversity As Revealed by the Population Structure Analysis of Germplasm Bank of Theobroma cacao L.

Beans of the species Theobroma cacao L., also known as cacao, are the raw material to produce chocolate. Colombian cacao has been classified as a fine flavor cacao that represents the 5% of cacao world's production. Colombian genetic resources from this species are conserved in ex situ and in-field germplasm banks, since T. cacao has recalcitrant seeds to desication and long-term storage. Currently, the collection of T. cacao of the Colombian Corporation of Agricultural Research (CORPOICA) has approximately 700 germplasm accessions. We conducted a molecular analysis of Corpoica's cacao collection and a morphological characterization of some accessions with the goal to study its genetic diversity and population structure and, to select interesting accessions for the cacao's breeding program. Phenotypic evaluation was performed based on 18 morphological traits and 4 biochemical traits. PCA analysis of morphological traits explained 60.6% of the total variation in seven components and 100% of the total variation of biochemical traits in four components, grouping the collection in 4 clusters for both variables. We explored 565 accessions from Corpoica's germplasm and 252 accessions from reference populations using 96 single nucleotide polymorphism (SNP) molecular markers. Molecular patterns of cacao Corpoica's collection were obtained amplifying specific alleles in a Fluidigm platform that used integrated circuits of fluids. Corpoica's collection showed highest genetic diversity [Expected Heterozygosity (HE = 0.314), Observed Heterozygosity (HO = 0.353)] that is reduced when reference populations were included in the dataset (HE = 0.294, HO = 0.261). The collection was divided into four clusters based on population structure analysis. Cacao accessions from distinct groups showed some taxonomic concordance and reflected their geographic origins. For instance, accessions classified as Criollo were clearly differentiated in one group and we identified two new Colombian genetic groups. Using a number of allelic variations based on 87 SNP markers and 22 different morphological/biochemical traits, a core collection with a total of 232 accessions was selected as a primary genetic resource for cacao breeders.

Association analysis for disease resistance to Fusarium oxysporum in cape gooseberry (Physalis peruviana L).

BACKGROUND: Vascular wilt caused by Fusarium oxysporum is the most important disease in cape gooseberry (Physalis peruviana L.) in Colombia. The development of resistant cultivars is considered one of the most cost-effective means to reduce the impact of this disease. In order to do so, it is necessary to provide breeders with molecular markers and promising germplasm for introgression of different resistance loci as part of breeding schemes. Here we described an association mapping study in cape gooseberry with the goal to: (i) select promising materials for use in plant breeding and (ii) identify SNPs associated with the cape gooseberry resistance response to the F. oxysporum pathogen under greenhouse conditions, as potential markers for cape gooseberry breeding. RESULTS: We found a total of 21 accessions with different resistance responses within a diversity panel of 100 cape gooseberry accessions. A total of 60,663 SNPs were also identified within the same panel by means of GBS (Genotyping By Sequencing). Model-based population structure and neighbor-joining analyses showed three populations comprising the cape gooseberry panel. After correction for population structure and kinship, we identified SNPs markers associated with the resistance response against F. oxysporum. The identification of markers was based on common tags using the reference genomes of tomato and potato as well as the root/stem transcriptome of cape gooseberry. By comparing their location with the tomato genome, 16 SNPs were found in genes involved in defense/resistance response to pathogens, likewise when compared with the genome of potato, 12 markers were related. CONCLUSIONS: The work presented herein provides the first association mapping study in cape gooseberry showing both the identification of promising accessions with resistance response phenotypes and the identification of a set of SNP markers mapped to defense/resistance response genes of reference genomes. Thus, the work also provides new knowledge on candidate genes involved in the P. peruviana - F. oxysporum pathosystem as a foundation for further validation in marker-assisted selection. The results have important implications for conservation and breeding strategies in cape gooseberry.

Genetic diversity and population structure in Physalis peruviana and related taxa based on InDels and SNPs derived from COSII and IRG markers.

The genus Physalis is common in the Americas and includes several economically important species, among them Physalis peruviana that produces appetizing edible fruits. We studied the genetic diversity and population structure of P. peruviana and characterized 47 accessions of this species along with 13 accessions of related taxa consisting of 222 individuals from the Colombian Corporation of Agricultural Research (CORPOICA) germplasm collection, using Conserved Orthologous Sequences (COSII) and Immunity Related Genes (IRGs). In addition, 642 Single Nucleotide Polymorphism (SNPs) markers were identified and used for the genetic diversity analysis. A total of 121 alleles were detected in 24 InDels loci ranging from 2 to 9 alleles per locus, with an average of 5.04 alleles per locus. The average number of alleles in the SNP markers was two. The observed heterozygosity for P. peruviana with InDel and SNP markers was higher (0.48 and 0.59) than the expected heterozygosity (0.30 and 0.41). Interestingly, the observed heterozygosity in related taxa (0.4 and 0.12) was lower than the expected heterozygosity (0.59 and 0.25). The coefficient of population differentiation FST was 0.143 (InDels) and 0.038 (SNPs), showing a relatively low level of genetic differentiation among P. peruviana and related taxa. Higher levels of genetic variation were instead observed within populations based on the AMOVA analysis. Population structure analysis supported the presence of two main groups and PCA analysis based on SNP markers revealed two distinct clusters in the P. peruviana accessions corresponding to their state of cultivation. In this study, we identified molecular markers useful to detect genetic variation in Physalis germplasm for assisting conservation and crossbreeding strategies.